function output = Qf2(lambda, p0, gamma)


output = lambda + 4*9.8*abs(p0)/(sqrt(lambda) + sqrt(lambda + 2*p0*gamma))

end
